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Summary. This article contains a theoretical overview of the physical properties 
of antiferromagnetic Mott insulators in spatial dimensions greater than one. Many 
such materials have been experimentally studied in the past decade and a half, 
and we make contact with these studies. The simplest class of Mott insulators have 
an even number oi S — 1/2 spins per unit cell, and these can be described with 
quantitative accuracy by the bond operator method: we discuss their spin gap and 
magnetically ordered states, and the transitions between them driven by pressure or 
an applied magnetic field. The case of an odd number of S = 1/2 spins per unit cell 
is more subtle: here the spin gap state can spontaneously develop bond order (so 
the ground state again has an even number of S* = 1/2 spins per unit cell), and/or 
acquire topological order and fractionalized excitations. We describe the conditions 
under which such spin gap states can form, and survey recent theories (T. Senthil 
et al, cond-mat/0312617) of the quantum phase transitions among these states and 
magnetically ordered states. We describe the breakdown of the Landau- Ginzburg- 
Wilson paradigm at these quantum critical points, accompanied by the appearance 
of emergent gauge excitations. 



1 Introduction 

The physics of Mott insulators in two and higher dimensions has enjoyed much 
attention since the discovery of cuprate superconductors. While a quantitative 
synthesis of theory and experiment in the superconducting materials remains 
elusive, much progress has been made in describing a number of antiferro- 
magnetic Mott insulators. A number of such insulators have been studied 
extensively in the past decade, with a few prominent examples being CaV409 
n, (C5Hi2N2)2Cu2Cl4 HElia, SrCu2(B03)2 TlCuCla 1313131101, 

and CS2CUCI4 |ll[ll2j . In some cases, it has even been possible to tune these 
insulators across quantum phase transitions by applied pressure |H] or by an 
apphed magnetic field |210I|7||nj- A useful survey of some of these experiments 
may be found in the recent article by Matsumoto et al. |l(Jj . 
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It would clearly be valuable to understand the structure of the global phase 
diagram of antiferromagnetic Mott insulators above one dimension. The com- 
pounds mentioned above would then correspond to distinct points in this 
phase diagram, and placing them in this manner should help us better under- 
stand the relationship between different materials. One could also classify the 
quantum critical points accessed by the pressure or field-tuning experiments. 
The purpose of this article is to review recent theoretical work towards achiev- 
ing this goal. We will focus mainly on the case of two spatial dimensions (d), 
but our methods and results often have simple generalizations to d = 3. 

One useful vantage point for opening this discussion is the family of Mott 
insulators with a gap to all spin excitations. All spin gap compounds dis- 
covered to date have the important property of being "dimerized" , or more 
precisely, they have an even number of 5 — 1/2 spins per unit cell In such 
cases, the spin gap can be understood by adiabatic continuation from the sim- 
ple limiting case in which the spins form local spin singlets within each unit 
cell. A simple approach that can be used for a theoretical description of such 
insulators is the method of bond operators |14l ITK) . This method has been 
widely applied, and in some cases provides an accurate quantitative descrip- 
tion of numerical studies and experiments |lfil ll()| . We will describe it here 
in Section 121 in the very simple context of a coupled dimer antiferromagnet; 
similar results are obtained in more complicated, and realistic, lattice struc- 
tures. Section[21will also describe the quantum phase transition(s) accessed by 
varying coupling constants in the Hamiltonian while maintaining spin rota- 
tion invariance (this corresponds to experiments in applied pressure): the spin 
gap closes at a quantum critical point beyond which there is magnetic order. 
Section O will discuss some of the important experimental consequences of 
this quantum criticality at finite temperatures. A distinct quantum critical 
point, belonging to a different universality class, is obtained when the spin 
gap is closed by an applied magnetic field — this is described in Section |3 

The remaining sections discuss the theoretically much more interesting and 
subtle cases of materials with an odd number of 5 = 1/2 spins per unit cell, 
such as La2Cu04 and CS2CUCI4. A complementary, but compatible, perspec- 
tive on the physics of such antiferromagnets may be found in the review article 
by Misguich and Lhuillier ^Jj- Antiferromagnets in this class can develop a 
spin gap by spontaneously breaking the lattice symmetry so that the lattice 
is effectively dimerized (see discussion in the following paragraph). There are 
no known materials with a spin gap in which the lattice symmetry has not 
been broken, but there is a theoretical consensus that spin gap states without 
lattice symmetry breaking are indeed possible in d > 1 ^Hl- The study of 
antiferromagnets with an odd number of 5* = 1/2 spins per unit cell is also 
important for the physics of the doped cuprates. These materials exhibit spin- 
gap-like behavior at low dopings, and many theories associate aspects of its 
physics with the spin gap state proximate to the magnetically ordered state 
of the square lattice antiferromagnet found in La2Cu04. 
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Section^will describe the nature of a spin gap state on the square lattice. 
We begin with the nearest-neighbor S — 1/2 Heisenberg Hamiltonian on the 
square lattice — this is known to have a magnetic Neel order which breaks 
spin rotation invariance. Now add further neighbor exchange couplings until 
magnetic order is lost and a spin gap appears. We will show that the ground 
state undergoes a novel, second-order quantum phase transition to a state with 
bond order: translational symmetry is spontaneously broken |19ll2()j so that the 
resulting lattice structure has an even number of S* = 1/2 spins per unit cell. 
So aspects of the non-zero spin excitations in this paramagnet are very similar 
to the "dimerized" systems considered in Section |21 and experimentally they 
will appear to be almost identical. Indeed, it may well be that the experimental 
materials initially placed in the class of Sectional are secretely systems in the 
class of Section 01 which have developed bond order driven by the physics of 
antiferromagets (as in Section 14. 1|) at some intermediate energy scale. The 
host lattice then distorts sympathetically to the bond order, and is effectively 
dimerized. Such materials will possess many more low-lying singlet excitations 
than those in the theory of Section |2 these excitations play an important 
role in the restoration of translational symmetry as we move towards the 
Neel state. Unfortunately, such singlet excitations are rather difhcult to detect 
experimentally. 

Section [S] will address the same issue as Section 0] but for the case of 
the triangular lattice. Here the spins are ordered in a non-collinear configu- 
ration in the magnetically ordered state, as is observed at low temperatures 
in CS2CUCI4 fJEl- We will argue that in this case there is a route to de- 
struction of magnetic order in which the resulting spin gap state preserves 
full lattice symmetry |^ |21]. Such a spin gap state has a novel 'topologi- 
cal' order which endows its excitations with charges under an emergent 
gauge force. Recent experimental measurements of the dynamic structure fac- 
tor of CS2CUCI4 appear to be described rather well by the excitations of this 
topologically ordered state at energies above which the magnetic order of the 
ground state emerges ^1 1^ . 

2 Coupled dimer antiferromagnet 

We begin by describing the quantum phase transition in a simple two- 
dimensional model of antiferromagnetically coupled S = 1/2 Heisenberg spins 
which has 2 spins per unit cell. The transition is tuned by varying a dimen- 
sionless parameter A. As we noted in Section ^ different 'dimerized' Mott 
insulators will correspond to different values of A, and the value of A can be 
tuned by applying pressure piTH]. 

We consider the "coupled dimer" Hamiltonian |2S] 



Hfj — J ^ ^ Si ■ Sj + XJ ^ ^ S,; • Sj 



(1) 
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Fig. 1. The coupled dimer antiferromagnet. Spins (S — 1/2) are placed on the 
sites, the A links are shown as full lines, and the B links as dashed lines. 

where Sj are spin- 1/2 operators on the sites of the coupled-ladder lattice 
shown in FigQ] with the A links forming decoupled dimers while the B links 
couple the dimers as shown. The ground state of Hd depends only on the 
dimensionless coupling A, and we will describe the low temperature (T) prop- 
erties as a function of A. We will restrict our attention to J > and < A < 1. 

Note that exactly at A = 1, Hd is identical to the square lattice antifer- 
romagnet, and this is the only point at which the Hamiltonian has only one 
spin per unit cell. At all other values of A Hd has a pair of S = 1/2 spins 
in each unit cell of the lattice. As will become clear from our discussion, this 
is a key characteristic which permits a simple theory for the quantum phase 
transition exhibited by Hd- Models with only a single S = 1/2 spin per unit 
cell usually display far more complicated behavior, and will be discussed in 
Sections 

We will begin with a physical discussion of the phases and excitations 
of the coupled dimer antiferromagnet, Hd in Section 12.11 We will propose 
a quantum field-theoretical description of this model in Section 1??^ we will 
verify that the limiting regimes of the field theory contain excitations whose 
quantum numbers are in accord with the phases discussed in Section [2. II and 
will then use the field theory to describe the quantum critical behavior both 
at zero and finite temperatures. 

2.1 Phases and their excitations 

Let us first consider the case where A is close to 1. Exactly at A = 1, Hd is 
identical to the square lattice Heisenberg antiferromagnet, and this is known to 
have long-range, magnetic Neel order in its ground state i.e. the spin-rotation 
symmetry is broken and the spins have a non-zero, staggered, expectation 
value in the ground state with 

(S,) - rj.Non, (2) 

where n is some fixed unit vector in spin space, rjj is ±1 on the two sublat- 
tices, and Nq is the Neel order parameter. This long-range order is expected 
to be preserved for a finite range of A close to 1. The low- lying excitations 
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Fig. 2. Schematic of the quantum paramagnet ground state for small A. The ovals 
represent singlet valence bond pairs. 
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Fig. 3. (a) Cartoon picture of the bosonic S = 1 excitation of the paramagnet. (b) 
Fission of the 5 = 1 excitation into two 5=1/2 spinous. The spinous are connected 
by a "string" of valence bonds (denoted by dashed ovals) which lie on weaker bonds; 
this string costs a finite energy per unit length and leads to the confinement of 
spinous. 



above the ground state consist of slow spatial deformations in the orienta- 
tion n: these are the familiar spin waves, and they can carry arbitrarily low 
energy i.e. the phase is 'gapless'. The spectrum of the spin waves can be ob- 
tained from a text-book analysis of small fluctuations about the ordered Neel 
state using the Holstein-Primakoff method such an analysis yields two 
polarizations of spin waves at each wavevector k = {kr^^ky) (measured from 
the antiferromagnetic ordering wavevector), and they have excitation energy 
= (c^fc^ + Cyfc^)"'^^^, with Cx,Cy the spin-wave velocities in the two spatial 
directions. 

Let us turn now to the vicinity of A = 0. Exactly at A = 0, Hd is the 
Hamiltonian of a set of decoupled dimers, with the simple exact ground state 
wavefunction shown in Fig [21 the spins in each dimer pair into valence bond 
singlets, leading to a paramagnetic state which preserves spin rotation in- 
variance and all lattice symmetries. Excitations are now formed by breaking 
a valence bond, which leads to a three-iold degenerate state with total spin 
5 = 1, as shown in Fig|2K- At A = 0, this broken bond is localized, but at fi- 
nite A it can hop from site-to-site, leading to a triplet quasiparticle excitation. 
Note that this quasiparticle is not a spin-wave (or equivalently, a 'magnon') 
but is more properly referred to as a spin 1 exciton or a triplon j27|. We pa- 
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Fig. 4. Ground states of Hd as a function of A Tiie quantum critical point is at 
Ac = 0.52337(3). Tlie compound TlCuCla undergoes a similar quantum phase 
transition under applied pressure 

rameterize its energy at small wavevectors k (measured from the minimum of 
the spectrum in the Brillouin zone) by 

^2 I „2 j,2 
'^-"^^ ^ ' 

where A is the spin gap, and Cx, Cy are velocities; we will provide an explicit 
derivation of © in Section Fig |21 also presents a simple argument which 
shows that the S = 1 exciton cannot fission into two S* = 1/2 'spinous'. 

The very distinct symmetry signatures of the ground states and excita- 
tions between A ~ 1 and A « make it clear that the two limits cannot 
be continuously connected. It is known that there is an intermediate second- 
order phase transition at [^HIEHI ^ — K = 0.52337(3) between these states as 
shown in Fig 21 Both the spin gap A and the Neel order parameter Nq vanish 
continuously as Ac is approached from cither side. 

2.2 Bond operators and quantum field theory 

In this section we will develop a continuum description of the low energy 
excitations in the vicinity of the critical point postulated above. There are a 
number of ways to obtain the same final theory: here we will use the method 
of bond operators |14l I15j , which has the advantage of making the connection 
to the lattice degrees of freedom most direct. We rewrite the Hamiltonian 
using bosonic operators which reside on the centers of the A links so that it 
is explicitly diagonal at A = 0. There are 4 states on each A link (|tT): ITi); 
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|J,t)j and and we associate these with the canonical singlet boson s and 
the canonical triplet bosons {a = x, y, z) so that 

|s)^st|o) = -^(|n)-UT)) ; K.)^4|0) = ^(|TT)-Ui)) ; 
|t,)^4|0) = -^(|TT> + Ui)) ; M ^4|0) = -i=(|n) + UT)). (4) 

Here |0) is some reference vacuum state which does not correspond to a phys- 
ical state of the spin system. The physical states always have a single bond 
boson and so satisfy the constraint 

sU + tlt^ = l (5) 

By considering the various matrix elements (s|Si|ic(), (s|S2|tQ), ■ • •, of the 
spin operators Si. 2 on the ends of the link, it follows that the action of Si and 
S2 on the singlet and triplet states is equivalent to the operator identities 

Slot — 2 ^■S^^Q + i^afi'ft^fjt^ , 

S2a = ^ {^-sha - 4s " i^al3-ft\t-^ ■ (6) 

where a,/3,7 take the values x,y,z, repeated indices are summed over and e is 
the totally antisymmetric tensor. Inserting jnjl into (^1, and using we find 
the following Hamiltonian for the bond bosons: 

Hd = Ho + Hi 



a{£, to) (tljrnaslnSe + ^a^maSmS^ + H.C.^ + b{i, to) 



[ieap-^tl^J^ll^ti^Sm + H.C.j + c{t,m) (tl^tl^aim/stip ~ tl^tl^ptjnatefi 



(7) 



where £, m label links in A, and a, h, c are numbers associated with the lattice 
couplings which we will not write out explicitly. Note that _ffi = at A = 0, 
and so the spectrum of the paramagnetic state is fully and exactly deter- 
mined. The main advantage of the present approach is that application of the 
standard methods of many body theory to , while imposing the constraint 
©, gives a very satisfactory description of the phases with A 7^ 0, including 
across the transition to the Neel state. In particular, an important feature of 
the bond operator approach is that the simplest mean field theory already 
yields ground states and excitations with the correct quantum numbers; so a 
strong fluctuation analysis is not needed to capture the proper physics. 
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A complete numerical analysis of the properties of in a self-consistent 
Hartree-Fock treatment of the four boson terms in Hi has been presented in 
Ref. |14j . In all phases the s boson is well condensed at zero momentum, and 
the important physics can be easily understood by examining the structure of 
the low energy action for the ta bosons. For the particular Hamiltonian 
the spectrum of the ta bosons has a minimum at the momentum (0, tt), and 
for large enough A the ta condense at this wavevector: the representation © 
shows that this condensed state is the expected Neel state, with the magnetic 
moment oscillating as in Q . The condensation transition of the ta is therefore 
the quantum phase transition between the paramagnetic and Neel phases of 
the coupled dimer antiferromagnet. In the vicinity of this critical point, we 
can expand the ta bose field in gradients away from the (0, tt) wavevector: so 
we parameterize 

i£,a(T) =i„(r,,r)e'^(0'^>'-^ (8) 

where t is imaginary time, r = (x, y) is a continuum spatial co-ordinate, and 
expand the effective action in spatial gradients. In this manner we obtain 



St = / d^rdT 



2 



+ - f (to^ta + H.C.) + Ki,\dJa\^ + Kiy\dyt. 

+ ^ {K2rid^taf + K2y{dytaf + H.C.) + • 



(9) 



Here C, Z?, Ki,2x.y are constants that are determined by the solution of the 
self-consistent equations, and the ellipses represent terms quartic in the ta- 
The action St can be easily diagonalized, and we obtain a 5' = 1 quasiparticle 
excitation with the spectrum 



£k 



iC + Kixkl + Kiyklf -{D + K2xkl + K2yklf] . (10) 



This is, of course, the triplon (or spin exciton) excitation of the paramagnetic 
phase postulated earlier in lO ; the latter result is obtained by expanding H10|l 
in momenta, with A — VC"^ — D"^ . This value of A shows that the ground 
state is paramagnetic as long as C > D, and the quantum critical point to 
the Neel state is at C = I?. 

The critical point and the Neel state are more conveniently described by 
an alternative formulation of St (although an analysis using bond operators 
directly is also possible [211) ■ It is useful to decompose the complex field ta 
into its real and imaginary parts as follows 

ta^ Z{ifa+iTTa), (H) 

where Z \s a. normalization chosen below. Insertion of Hll|l into Q shows that 
the field tTq, has a quadratic term ~ (C -I- D)ti^, and so the co-efficient of tt^ 
remains large even as the spin gap A becomes small. Consequently, we can 
safely integrate tt^ out, and the resulting action for Lpa takes the form 
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Sip = j d rdr 



(12) 

Here we have chosen Z to fix the co-efficient of the temporal gradient term, 
and s^C^ - D^. 

The formulation makes it simple to explore the physics in the region 
s < 0. It is clear that the effective potential of fa has a minimum at a non- 
zero (^Q, and that {i-Pa) oc iVo, the Neel order parameter in (j^J. It is simple to 
carry out a small fluctuation analysis about this saddle point, and we obtain 
the doublet of gapless spin-wave modes advertised earlier. 

We close this subsection by noting that all of the above results have a direct 
generalization to other lattices, and also to spin systems in three dimensions. 
Matsumoto et al. ^U] have applied the bond operator method to TlCuCla 
and obtained good agreement with experimental observations. One important 
difference that emerges in such calculations on some frustrated lattices [30] 
is worth noting explicitly here: the minimum of the spectrum need not 
be at special wavevector like (0, tt), but can be at a more generic wavevector 
Q such that Q and —Q are not separated by a reciprocal lattice vector. A 
simple example which we consider here is an extension of ^ in which there 
are additional exchange interactions along all diagonal bonds oriented 'north- 
east' (so that the lattice has the connectivity of a triangular lattice). In such 
cases, the structure of the low energy action is different, as is the nature of 
the magnetically ordered state. The parameterization JHJ must be replaced by 

t^T) = hain, r)e^Q"-^ + haire. r)e-^Q-'-^ (13) 

where ^1^20 are independent complex fields. Proceeding as above, we find that 
the low energy effective action (|12() is replaced by 



Sd, — d rdr 



2 

2 



1(1 



(14) 



where now is a complex field such that ~ {tia) ~ (4a)- Notice that 
there is now a second quartic term with co-efficient v.li v > 0, configurations 
with = are preferred: in such configurations = ^T-ia + *^2a, where 
'^i,2q are two equal-length orthogonal vectors. Then from l|13() and ^ it is 
easy to see that the physical spins possess spiral order in the magnetically 
ordered state in which (p^ is condensed. A spiral state is illustrated in Fig llHI 
and we will have more to say about this state in Sectional For the case w < 0, 
the optimum configuration has 'Pa = UaC^^ where Ua is a real vector: this 
leads to a magnetically ordered state with spins polarized collinearly in a spin 
density wave at the wavevector Q. 



10 Subir Sachdev 



2.3 Quantum criticality 

We will restrict our discussion here to the critical point described by 5^. 
Similar results apply to S$ for the parameter regime in which it exhibits a 
second order transition |31| . ExperimentaUy, the results below are relevant 
to materials that can be tuned across the magnetic ordering transition by 
applied pressure (such as TlCuCla [H]), or to materials which happen to be 
near a critical point at ambient pressure (such as LaCu02.5 '32' ) . 

The field theory 5^ is actually a familiar and well-studied model in the 
context of classical critical phenomena. Upon interpreting r as a third spatial 
co-ordinate, becomes the theory of a classical 0(3)-invariant Heisenberg 
ferromagnet at finite temperatures (in general a d dimensional quantum anti- 
ferromagnet will map to a d + 1 dimensional classical Heisenberg ferromagnet 
at finite temperature The Curie transition of the Heisenberg ferromag- 

net then maps onto the quantum critical point between the paramagnetic 
and Neel states described above. A number of important implications for the 
quantum problem can now be drawn immediately. 

The theory 5^ has a 'relativistic' invariance, and consequently the dynamic 
critical exponent must be z = 1. The spin correlation length will diverge at the 
quantum critical point with the exponent v = 0.7048(30). The spin gap 
of the paramagnet, A, vanishes as Z\ ~ {Xc — X)^'^ , and this prediction is in ex- 
cellent agreement with the numerical study of the dimerized antiferromagnet 



A somewhat more non-trivial consequence of this mapping is in the struc- 
ture of the spectrum at the critical point A = Ac. At the Curie transition of the 
classical ferromagnet it is known 35 that spin correlations decay as ~ l/p'^~'^, 
where p is the 3-component momentum in the 3-dimensional classical space. 
We can now analytically continue this expression from its Pz dependence in 
the third classical dimension to the real frequency, uj, describing the quan- 
tum antiferromagnet. This yields the following fundamental result for the 
dynamic spin susceptibility, x(/c,w), at the T ^ quantum critical point of 
the coupled-dimer antiferromagnet: 



where e is a positive infinitesimal. Note that in (|15|l the momentum k is 
measured from the (tt, tt) ordering wavevector of the Neel state. The exponent 
77 is the same as that of the classical Heisenberg ferromagnet, and has a rather 
small value |34| : rj « 0.03. However, the non-zero r/ does make a significant 
difference to the physical interpretation of the excitations at the critical point. 
In particular note that lmx{k,uj) does not have a pole at any fc, but rather a 
continuum spectral weight above a threshold energy 1^ 



m 



1 



(15) 





/„2 1,2 I „2 1,2 \ 

.2 1,2 _ „2 l.2\ 1"''/' 
■x "'X ^y '^y J 



(16) 
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where 9 is the unit step function. This indicates there are no quasiparticles at 
the critical point, and only a dissipative critical continuum. 

There is also some very interesting structure in the quantum critical dy- 
namic response at nonzero T [56, 37 . Here, one way to understand the physics 
is to approach the critical point from the paramagnetic side (A < Ac). As we 
noted earlier, the paramagnetic phase has well-defined 'triplon' or 'spin exci- 
ton' excitations ta, and these have an infinite lifetime at T = 0. At T > 0, ther- 
mally excited ta quasiparticles will collide with each other via their scattering 
amplitude, m, and this will lead to a finite lifetime |^|nHl- Now approach 
A = Ac- The renormalization group analysis of tells us that the quartic 
coupling u approaches a fixed point value in the critical region. This means 
that u is no longer an arbitrary parameter, and an appropriately defined ta 
scattering amplitude must also acquire universal behavior. In particular, the 
ta lifetime is determined by the only energy scale available, which is fcsT. So 
we have the remarkable result that the characteristic spin relaxation time is 
a universal number times h/{kBT). More precisely, we can write for the local 
dynamic spin susceptibility xl{^) — J (Pkx{k,u>) the universal scaling form 

Imx.M=r'F(^) (17) 

Here F is a universal function which has the limiting behaviors 

^^^)^{^gn(c.)M<|fe|5l • (1^) 

Note that F has a smooth linear behavior in the regime \huj\ ^ fcsT, and 
this is similar to any simple dissipative system. The difference here is that the 
co-efficient of dissipation is determined by ksT alone. 

The quantum critical behavior described here is expected to apply more 
generally to other correlated electron systems, provided the critical theory has 
non-linear couplings which approach fixed point values. 

3 Influence of an applied magnetic field 

An important perturbation that can be easily applied to antiferromagnets in 
the class discussed in Section|21is a uniform magnetic field. The Zeeman energy 
in available fields can often be comparable to the typical antiferromagnetic 
exchange constant J, and so the ground state can be perturbed significantly. 
It is therefore of interest to understand the evolution of the phase diagram in 
Fig 0] under an applied field of arbitrary strength. 

We are interested here in the evolution of the ground state as a function 
of B where the Hamiltonian Hd in Q is transformed as 



(19) 
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Most of the basic features can actually be understood quite easily in a simple 
extension of the self-consistent Hartree-Fock theory of bond bosons that was 
discussed in Section Under the transformation (|19|l . it is easily seen from 
® that 

Hd ^ Hd + iBa ^ Cafljtlpte^ (20) 

eeA 

The presence of a non-zero B breaks spin rotation invariance and so all the 
self-consistent expectation values of operator bilinears have to reflect this 
reduced symmetry in the Hartree-Fock theory. Apart from this the mechanics 
of the computation mostly remain the same. However, for stronger fields, 
it is sometimes necessary to allow for broken translational symmetry in the 
expectation values, as the ground state can acquire a modulated structure. 

We will discuss the results of such an analysis in weak and strong fields in 
the following subsections. 

3.1 Weak fields 

For weak fields applied to the paramagnet (specifically, for fields B < A, the 
influence of (|2()|l can be understood exactly. The coupling to B involves an 
operator which commutes with the remaining Hamiltonian (the total spin), 
and hence the wavefunction of the ground state remains insensitive to the 
value of B. The same applies to the wavefunctions of the excited states. How- 
ever, the excited states can have non-zero total spin and so their energies do 
depend upon B. In particular the triplet ta quasiparticle with energy or 
pUI) carries total spin S — 1, and consequently we conclude that this triplet 
splits according to 

£k £fe - rriB (21) 

with m = 0, ±1. Note that the lowest energy quasiparticle (with m = 1) has a 
positive energy as long as B < A, and this is required for the stability of the 
paramagnet. So the phase boundary of the paramagnetic phase is exactly B — 
A, and using A ~ (Ac — A)^'', we can sketch the boundary of the paramagnetic 
phase as in Fig El 

What happens beyond the paramagnetic phase? As in Section 12.21 we 
answer this question by using the transformation (|11() , and by examining the 
analog of S^p under a non-zero B. Using H20(l . and integrating out tTq,, we now 
find that the action S^p in (|12(l remains unchanged apart from the mapping 

na 

{dr^Pa) ^ {dr(Pa+iiaf3-fBf3ip~f) . (22) 

The action (O, ^ can now be analyzed by a traditional small fluctuation 
analysis about ipa — 0. Let us assume that B = (0, 0, B) is oriented along the 
z axis. Then the co-efficient of (pi is s, while that of (p^. + ipy is s — B^. This 
suggests that we focus only on the components of ipa in the plane orthogonal 
to B, and integrate out the component of ipa along the direction of B. Indeed, 
if we define 
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;^ 1 

Fig. 5. Evolution of the phases of Fig |1| under a weak field B (magnetization 
plateau at large B, appearing in Fig|S] are not shown). The paramagnetic phase has 
exactly the same ground state wavefunction as that at B = 0. The phase boundary 
behaves like B ~ (Ac — A)^". The B field is oriented vertically upwards, and the static 
moments in the canted phase can rotate uniformly about the vertical axis. The phase 
boundary at non-zero B is described by the 2 = 2 dilute Bose gas quantum critical 
theory. The phase diagram of TlCuCls in applied pressure and magnetic field looks 
similar to the one above |1U| . The corresponding phase diagram of the field-induced 
magnetic ordering transition of a superconductor (rather than a Mott insulator) has 
been investigated recently and successfully applied to experiments on the doped 
cuprates; this phase diagram of the superconductor has significant differences from 
the one above. 



and integrate out cpz, then we obtain from l|12(l . (|22|l the effective action for 



= I (frdr 



2 , S lo .7-,|2 . i.7-,i2 , " l.T,l4 

(24) 

Here, n = {s — B^)/2B, and we have retained only leading order temporal 
and spatial gradients and the leading dependence of u. Clearly, this is the 
theory of a Bose gas in the grand canonical ensemble at a chemical potential 
/K, with a repulsive short-range interaction 01]. At T = 0, and A* < 0, such a 
theory has a ground state which is simply the vacuum with no Bose particles. 
Here, this vacuum state corresponds to the spin gap antiferromagnet, and the 
^-independence of the ground state of the antiferromagnet corresponds here 
to the /X independence of the ground state of S,p. There is an onset of a finite 
density of bosons in S,p for > 0, and this onset therefore corresponds to the 
quantum phase transition in the antiferromagnet at -B = Z\. So we must have 
= in iSij/ at precisely the point where B = A: the value of /i quoted above 
shows that this is true at zeroth order in u, and higher order terms in u must 
conspire to maintain this result. 
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The above analysis makes it clear that the /i > region of Sxp will describe 
the quantum phase transition out of the paramagnet at non-zero B. This 
transition is merely the formation of a Bose- Einstein condensate of the m = 1 
component of the triplon bosons. For fi > we have a finite density of ^ 
bosons which Bose condense in the ground state, so that {^) ^ 0. From H23() 
we see that this Bose condensation corresponds to antiferromagnetic order in 
the plane perpendicular to B. Indeed, the phase of this Bose condensate is 
simply the orientation of the spins in the x, y plane, and so here this phase is 
directly observable. Further, by taking derivatives of (|19|) and w.r.t. B, we 
see that the density of bosons is proportional to the magnetization per spin, 
/2, in the direction parallel to B: 

^^^Y.(s,.)^m'), (25) 

i 

where N is the total number of spins. Consequently, the average magnetic 
moments in the non-paramagnetic phase are in a 'canted' configuration, as 
shown in FigEl The quantum phase transition between the paramagnet and 
the canted state is described by the theory of the density onset in a Bose gas: 
this theory has z = 2, v = 1/2, and an upper critical dimension of d = 2 

miiT]. 

We conclude this section by noting that interesting recent work 031 has 
examined the Bose-Einstein condensation of the m — 1 triplon bosons in a 
random potential. This is achieved by studying Tli-^^KajCuCla, where the sto- 
ichiometric disorder among the non-magnetic ions acts as a random potential 
on the triplons. 



3.2 Strong fields 

We have seen above that applying a magnetic field eventually leads to the 
onset of a ferromagnetic moment in the directions of the applied field. How 
does this moment evolve as we continue to increase the field? Eventually, B 
will become so large that it pays to have all the spins polarized in the direction 
of the field: this corresponds to a saturation in the magnetization, and making 
B even stronger will not change the ground state. In terms of the t bosons, 
this fully polarized state, \FP), with /2 = 1/2, is seen from (12011 or Q) to 
correspond exactly to 

\FP)^l[^^k^\0). (26) 

So there must be at least one more quantum phase transition as a i? is in- 
creased: this is transition from the \FP) state at very large _B to a state with 
a continuously varying ferromagnetic moment which eventually reaches the 
saturation value from below. 
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A theory for the transition away from the \FP) state with decreasing B 
can be developed using methods very similar to those used in Section 12.21 
and 13. II We treat the quartic terms in Q) in a Hartree-Fock approximation, 
and examine small fluctuations away from the \FP) state. These are domi- 
nated by excitation which create tz quanta (which have m = 0) on the dimers, 
and so the effective theory is expressed in terms of 

^tl{t,~zty). (27) 

Indeed, it is not difficult to see that the resulting theory for ^ has exactly 
the same form as (|24ll . Now the ^ for decreases with increasing B, and we 
have /i = at the critical field at which \FP) first becomes the ground state. 
Furthermore, now measures the deviation away from ^2 = 1/2. Apart 

from this 'inversion' in the field axis, it is clear that the universality class of 
the present transition is identical to that discussed in Section |0] 

A further possibility for a plateau in the value of f2 with increasing B 
is worth mentioning (Sj, as analogs are realized in SrCu2 (603)2 0Hj and 
NH4CUCI3 0ni- So far we have found plateaus at i7 = for i? < zi, and at i7 = 
1/2 for large B. For the J7 = 1/2 state we had every dimer with a {fl+if^y) /^/2 
boson. Now imagine that these bosons form a Wigner-crystalline state so that 
there are p such bosons for every q dimers; here < p < q, q > 1, are integers. 
Such a state will have f2 = p/{2q), and breaks the translational symmetry of 
the underlying dimer antiferromagnet such that there are q dimers per unit 
cell (or 2q spins per unit cell). The energy gap towards boson motion in the 
Wigner crystal {i.e. its incompressibility) will ensure that i? is stable under 
small variations of i?. In this manner we can obtain a magnetization plateau 
at f2 = p/{2q) in a state with a unit cell of q dimers. 

We summarize the considerations of this subsection in Fig (SJ showing a 
possible evolution of i7 in a model similar to Hd in As we have already 
noted, the plateau onset transitions at f2 = and fl = 1/2 are both described 
by the z = 2 dilute Bose gas theory H24|l . The transitions in and out of other 
fractional plateaus are potentially more complicated because these involve 
spontaneous breaking of translational symmetry. The translation symmetry 
could be restored at the same point at which there is onset of superfluid 
order — this is possibly a first order transition with a jump in the value of f2. 
Alternatively, there could be an intermediate 'supersolid' phase, in which case 
the plateau transition has the same broken translational symmetry on both 
sides of it, placing it also in the class of H24|l . 

4 Square lattice antiferromagnet 

This section will address the far more delicate problem of quantum phase 
transitions in antifcrromagnets with an odd number oi S = 1/2 spins per unit 
cell. We will mainly concern ourselves with square lattice Hamiltonians of the 
form 
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Fig. 6. Magnetization density, i7, defined in 1251 as a function of the applied mag- 
netic field. The plateau shown at = is present provided the zero field state is a 
paramagnet i.e. X < Ac. The full saturation plateau at J7 = 1/2 is always present. 
The plateau at i7 = 1/4 is not present in the nearest-neighbor model Ha in but 
it is believed that such a plateau will appear upon including frustrating exchange in- 
teractions; this plateau will involve a broken translational symmetry in the coupled 
dimer antiferromagnet. Such magnetization plateaux are found in SrCu2 (603)2 |45| 
and NH4CUCI3 gS] 



Here J is a nearest-neighbor antiferromagnetic exchange and the eUipses rep- 
resent further short-range exchange interactions (possibly involving multiple 
spin ring exchange) which preserve the full symmetry of the square lattice. 
The model is a member of the class Hg only at A = 1; at other values 
of A the symmetry group of the square lattice is exphcitly broken, and the 
doubhng of the unit cell was crucial in the analysis of Section El With full 
square lattice symmetry, the paramagnetic phase is not determined as simply 
as in the small A expansion, and we have to account more carefully for the 
'resonance' between different valence bond configurations. 

One ground state of Hg is, of course, the Neel state characterized by 
this is obtained in the absence of the interactions denoted by ellipses in (^5)1 . 
Now imagine tuning the further neighbor couplings in (|28|) so that spin rota- 
tion invariance is eventually restored and we obtain a paramagnetic ground 
state. We can divide the possibilities for this state into two broad classes, 
which we discuss in turn. 

In the first class of paramagnets, no symmetries of the Hamiltonian are 
broken, and the spins have paired with each other into valence bond singlets 
which strongly resonate between the large number of possible pairings: this is 
a resonating valence bond (RVB) liquid 071I3H1- We will discuss such states 
further in Section (SJ they have a connection with magnetically ordered states 
with non-coUinear magnetic order, unlike the collinear Neel state of the nearest 
neighbor square lattice antiferromagnet. 




(28) 
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In the second class of paramagnets, the valence bond singlets sponta- 
neously crystallize into some configuration which necessarily breaks a lat- 
tice symmetry. A simple example of such a bond- ordered paramagnet is the 
columnar state we have already considered in Fig |21 For the dimerized an- 
tiferromagnet Hd, the bond configuration in Fig [21 was chosen explicitly in 
the Hamiltonian by the manner in which we divided the links into classes A 
and B for A 7^ 1. For Hg, there is no such distinction between the links, and 
hence a state like Fig |21 spontaneously breaks a lattice symmetry. Further- 
more, there are 3 other equivalent states, obtained by successive 90 degree 
rotations of Fig[21 about any lattice site, which are completely equivalent. So 
for Hs, the bond-ordered paramagnet in Figl^jis four- fold degenerate. Going 
beyond simple variational wavefunctions like Fig El the bond-ordered states 
are characterized by a bond order parameter 



the values of on the links of the lattice in a bond-ordered state have a 
lower symmetry than the values of the exchange constants Jij in the Hamil- 
tonian. We will develop an effective model for quantum fluctuations about 
the coUinear Neel state in Hg below, and will find that such bond-ordered 
paramagnets emerge naturally 

Let us now try to set up a theory for quantum fluctuations about the Neel 
state (0). It is best to do this in a formulation that preserves spin rotation 
invariance at all stages, and this is facilitated by the coherent state path 
integral (see Chapter 13 of Ref. 01]). The essential structure of this path 
integral can be understood simply by looking at a single spin in a magnetic 
field h with the Hamiltonian Hi = — h • S. Then its partition function at a 
temperature T is given by 



Here S is the angular momentum of the spin S (we are interested primarily 
in the case S = 1/2) and n(T) is a unit 3-vector with n(0) = n(l/T). So the 
above path integral is over all closed curves on the surface of a sphere. The 
first term in the action of the path integral is the crucial Berry phase: ^[n(r)] 
is half the oriented area enclosed by the curve n(r) (the reason for the half 
will become clear momentarily). Note that this area is only defined modulo 
47r, the surface area of a unit sphere. The expression H3U|) has an obvious 
generalization to the lattice Hamiltonian Hg: the action adds up the Berry 
phases of every spin, and there is an additional energy term which is just the 
Hamiltonian with the replacement Sj Suj . 

We are now faced with the problem of keeping track of the areas enclosed 
by the curves traced out by all the spins. This seems rather daunting, partic- 
ularly because the half-area y^[n(r)] is a global object defined by the whole 
curve, and cannot be obviously be associated with local portions of the curve. 



Q 




(29) 
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Fig. 7. The path traced out by a single spin on the unit sphere in imaginary time. 
After discretizing time, the area enclose by the path is written as the sum over the 
areas of spherical triangles: Aij is half the area of the triangle with vertices no, n^, 
Uj . Different choices for the arbitrary point no correspond to different gauge choices 
associated with 13211 and 13411 . 

One convenient way to proceed is illustrated in Fig {7\ discretize imaginary 
time, choose a fixed arbitrary point no on the sphere, and thus write the area 
as the sum of a large number of spherical triangles. Note that each triangle is 
associated with a local portion of the curve n(T). 

We now need an expression for y^(ni, n2, na), defined as half the area of 
the spherical triangle with vertices ni, n2, n^. Complicated expressions for 
this appear in treatises on spherical trigonometry, but a far simpler expression 
is obtained after transforming to spinor variables [50j. Let us write 



where a,b =1,1 and we will always assume an implied summation over such 
indices, (Tab are the Pauli matrices, and Zj-^, Zj^ are complex numbers obeying 
p + \zji\^ = 1. Note that knowledge of Uj only defines zja up to a U(l) 
gauge transformation under which 



(31) 



z. 




(32) 



Then, associated with each pair of vertices n.^, tij we define 



A-ij — arg [z^^Zja 
Under the gauge transformation H32II we have 



(33) 



(34) 
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i.e. Aij behaves like a U(l) gauge field. Note also that Aij is only defined 
modulo 2iT, and that Aji = —Aij. For future use, we also mention the following 
identity, which follows from (|31|) and H33|) : 

1/2 

e*-^'^ . (35) 



The classical result for the half-area of the spherical triangle can be written 
in the simple form in terms of the present U(l) gauge variables: 

^(ni,n2,n3) =^i2+^23+-43i (36) 

We chose ^ as a half-a.rea. earlier mainly because then the expressions (|33f) 
and (|35|l come out without numerical factors. It is satisfying to observe that 
this total area is invariant under (|34|l . and that the half-area is ambiguous 
modulo 27r. 

Using H3fi|l. we can now write down a useful expression for ^[n(T)]. We 
assume that imaginary time is discretized into times tj separated by intervals 
At. Also, we denote by j + r the site at time tj + At, and define Ajj+r = Ajt ■ 
Then 

^[n(r)]=^^,. (37) 

3 

Note that this expression is a gauge-invariant function of the U(l) gauge field 
AjT, and is analogous to the quantity sometimes called the Polyakov loop. 

We are now ready to write down the first form proposed effective action 
for the quantum fluctuating Neel state. We do need to address some simple 
book-keeping considerations first: 

(i) Discretize spacetime into a cubic lattice of points j. Note that the same 
index j referred to points along imaginary time above, and to square lattice 
points in Hg. The meaning of the site index should be clear from the context. 
(m) On each spacetime point j, we represent quantum spin operator S^- by 

S, = r^.Sn,, (38) 

where iij is a unit vector, and rjj — ±1 is the sublattice staggering factor 
appearing in Q . This representation is that expected from the coherent state 
path integral, apart from the rjj factor. We have chosen to include rjj because of 
the expected local antiferromagnetic correlations of the spins. So in a quantum 
fluctuating Neel state, we can reasonably expect rij to be a slowly varying 
function of j. 

{in) Associated with each nj, define a spinor zja by (|^ . 
(iv) With each link of the cubic lattice, we use (|33|l to associate with it a 
Aj^ — Here fi = x,y,T extends over the 3 spacetime directions. 

With these preliminaries in hand, we can motivate the following effective 
action for fluctuations under the Hamiltonian H^: 
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J" 3 y (ij) 3 / 



.(39) 



Here the summation over {if) extends over nearest neighbors on the cubic 
lattice. The integrals are over the Zja-, and the rij and Ajr are dependent 
variables defined via (|31|l and H33() . Note that both terms in the action are 
invariant under the gauge transformation H32|l: consequently, we could equally 
well have rewritten Z as an integral over the , but it turns out to be more 
convenient to use the Zj^ and to integrate over the redundant gauge degree of 
freedom. The first term in the action contains the energy of the Hamiltonian 
Hs , and acts to prefer nearest neighbor rij which are parallel to each other — 
this "ferromagnetic" coupling between the rij in spacetime ensures, via 
that the local quantum spin configurations are as in the Neel state. The second 
term in the action is simply the Berry phase required in the coherent state 
path integral, as obtained from (|30(l and H37(l : the additional factor of 77^ com- 
pensates for that in (|38|l . The dimensionless coupling g controls the strength 
of the local antiferromagnetic correlations; it is like a "temperature" for the 
ferromagnet in spacetime. So for small g we expect Z to be in the Neel phase, 
while for large g we can expect a quantum- "disordered" paramagnet. For a 
much more careful derivation of the partition function Z from the underlying 
antiferromagnet iJg, including a quantitative estimate of the value of 5, see 
e.g. Chapter 13 of Ref. 

While it is possible to proceed with the remaining analysis of this section 
using Z, we find it more convenient to work with a very closely related alterna- 
tive model. Our proposed theory for the quantum fluctuating antiferromagnet 
in its final form is [511 1^ 




1 



2j+M,a + C.C.) -f i2S ^ ?7j Ajr j . (40) 

Note that we have introduced a new field Ajy^, on each link of the cubic lattice, 
which is integrated over. Like ^i^, this is also a U(l) gauge field because all 
terms in the action above are invariant under the analog of H34II : 

^ ^ "/"i + "^^i+A'- (41) 

The very close relationship between Z and Z may be seen |51j by explicitly 
integrating over the in l|4()|l : this integral can be done exactly because the 
integrand factorizes into terms on each link that depend only on a single Ajy^. 
After inserting (|35|l into (|40|l . the integral over the j/i link is 

r ^ exp ^2(1 + ".-".+/^))^/^ ^^3(^^.^ _ ^ ^2S^,^^.A,^ 
.In l-K \ q / 
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= hss^ 



(2(1 + n, • n,+,J) 



1/2 



exp(i2577jV^-^) (42) 



where the resuh involves either the modified Bessel function /q (for /i — x, y) 
or /2s (for n = t). We can use the identity 142|) to perform the integral over 
Ajfj, on each link of H40|l . and so obtain a partition function, denoted Z' , as 
an integral over the Zja only. This partition function Z' has essentially the 
same structure as Z in l|39|l . The Berry phase term in Z' is identical to that 
in Z. The integrand of Z' also contains a real action expressed solely as a 
sum over functions of • n^- on nearest neighbor links: in Z this function is 
simply Ui ■ Uj/g, but the corresponding function obtained from is more 
complicated (it involves the logarithm of a Bessel function), and has distinct 
forms on spatial and temporal links. We do not expect this detailed form of 
the real action function to be of particular importance for universal properties: 
the initial simple nearest-neighbor ferromagnetic coupling between the xij in 
H39() was chosen arbitrarily anyway. So we may safely work with the theory Z 
in (|4()|l henceforth. 

One of the important advantages of H4()|l is that we no longer have to keep 
track of the complicated non- linear constraints associated with 1)31(1 and ((33|l : 
this was one of the undesirable features of (|39|l . In Z, we simply have free 
integration over the independent variables Zja and Aj^ . The remainder of this 
section will be devoted to describing the properties of as a function of the 
coupling g. 

The theory Z in 1(40(1 has some resemblance to the so-called CP^^^ model 
from the particle physics literature [531 1541 1^ : our indices a,b take only 2 
possible values, but the general model is obtained when a,b — 1 . . .N, and 
we will also find it useful to consider Z for general N. The case of general 

describes SU(A^) and Sp(A^) antiferromagnets on the square lattice (19( . 
Note also that it is essential for our purposes that the theory is invariant 
under Aj^^ Aj^ + 27r, and so the U(l) gauge theory is compact. Finally 
our model contains a Berry phase term (which can be interpreted as a J^A^ 
term associated with a current Jj^^ — 2Srij5^r of static charges ±25* on each 
site) which is not present in any of the particle physics analyses. This Berry 
phase term will be an essential central actor in all of our results below for the 
paramagnetic phase and the quantum phase transition. 

The properties of Z are quite evident in the limit of small g. Here, the par- 
tition function is strongly dominated by configurations in which the real part 
of the action is a minimum. In a suitable gauge, these are the configurations in 
which Zja = constant, and by ((31(1 . we also have Uj a constant. This obviously 
corresponds to the Neel phase with (O . A Gaussian fluctuation analysis about 
such a constant saddle point is easily performed, and we obtain the expected 
spectrum of a doublet of gapless spin waves. 

The situation is much more complicated for large g where we should nat- 
urally expect a paramagnetic phase with (Sj) — (rij) = 0. This will be dis- 
cussed in some detail in Section 14.11 Finally, we will address the nature of 
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the quantum phase transition between the Neel and paramagnetic phases in 
Section lO 

4.1 Paramagnetic phase 

The discussion in this section has been adapted from another recent review 
by the author 

For large we can perform the analog of a 'high temperature' expansion 
of Z in H4Ufl . We expand the integrand in powers of 1/g and perform the 
integral over the Zja term-by-term. The result is then an effective theory for 
the compact U(l) gauge field Aj^ alone. An exphcit expression for the effective 
action of this theory can be obtained in powers oil/ g: this has the structure of 
a strong coupling expansion in lattice gauge theory, and higher powers of 1/g 
yield terms dependent upon gauge-invariant U(l) fluxes on loops of all sizes 
residing on the links of the cubic lattice. For our purposes, it is sufficient to 
retain only the simplest such term on elementary square plaquettes, yielding 
the partition function 



where e^i^A is the totally antisymmetric tensor in three spacetime dimensions. 
Here the cosine term represents the conventional Maxwell action for a compact 
U(l) gauge theory: it is the simplest local term consistent with the gauge 
symmetry (|41f) and which is periodic under —t Ajf^ + 27r; closely related 
terms appear under the l/g expansion. The sum over □ in (|43|) extends over 
all plaquettes of the cubic lattice, is the standard discrete lattice derivative 
{Af^fj = fj^f^ — fj for any /j), and is a couphng constant. We expect the 
value of e to increase monotonically with g. 

As is standard in duality mappings, we first rewrite the partition function 
in 2 -|- 1 spacetime dimensions by replacing the cosine interaction in 143|l by a 
Villain sum |5S1|S7] over periodic Gaussians: 



where the q-jp, are integers on the links of the dual cubic lattice, which pierce 
the plaquettes of the direct lattice. Throughout this article we will use the 
index j to refer to sites of this dual lattice, while j refers to the direct lattice 
on sites on which the spins are located. 

We will now perform a series of exact manipulations on (|44|l which will 
lead to a dual interface model |19[ I2(J[|^ . This dual model has only positive 




(43) 




(44) 
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Fig. 8. Specification of the non-zero values of the fixed field Oj^. The circles are the 
sites of the direct lattice, j, while the crosses are the sites of the dual lattice, j; the 
latter are also offset by half a lattice spacing in the direction out of the paper (the 
jj, = T direction). The a?^ are all zero for n = t,x, while the only non-zero values of 
Oijy are shown above. Notice that the a° flux obeys 1461 . 



weights — this fact, of course, makes it much more amenable to a standard 
statistical analysis. This first step in the duahty transformation is to rewrite 
H44fl by the Poisson summation formula: 



= XI ^'^P ( ~y m ~ * X! <^t^'^xo,jt^^'^^jX ) : (45) 
{a,^} \ 3 a J 

where aj^ (like qj^) is an integer- valued vector field on the hnks of the dual 
lattice (here, and below, we drop overall normalization factors in front of the 
partition function). Next, we write the Berry phase in a form more amenable to 
duality transformations. Choose a 'background' a^^ — a?^ flux which satisfies 

e^^xA^Qjx = %(5^r, (46) 

where j is the direct lattice site in the center of the plaquette defined by the 
curl on the left-hand-side. Any integer-valued solution of (|46l) is an acceptable 
choice for a?^ , and a convenient choice is shown in Fig|Sl Using H46|l to rewrite 
the Berry phase in 144(1 . applying H45|l . and shifting aj^ by the integer 25a?^, 
we obtain a new exact representation of Za in H44|) : 




^-=En/„ ^exp(-^^(a,,-25aP 



2 



□ 



^uA,x ■ (47) 
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The integral over the Aj^ can be performed independently on each link, and 
its only consequence is the imposition of the constraint e^j^i^xA^ajx — 0. We 
solve this constraint by writing aj^ as the gradient of a integer-valued 'height' 
hj on the sites of the dual lattice, and so obtain 

= ^ exp ( Y.^^^'^S - 2Sal)A ■ (48) 

We emphasize that, apart from an overall normalization, we have Zh — Za 
exactly. This is the promised 2+1 dimensional interface, or height, model in 
almost its final form. 

The physical properties of (|48|) become clearer by converting the "frustra- 
tion" in into offsets for the allowed height values. This is done by 
decomposing a?^ into curl and divergence free parts and writing it in terms 
of new fixed fields, Xj and 3^j^ as follows: 

a^"^ = A^Xj + e^^^xA^y^x. (49) 

The values of these new fields are shown in FigO Inserting H49|l into H48|> . we 
can now write the height model in its simplest form j20| 

5^exp(-^^(Z\,i//), (50) 

{H,} \ 3 J 

where 

= hj - 2SXj (51) 

is the new height variable we shall work with. Notice that the have 
dropped out, while the Xj act only as fractional offsets (for S not an even 
integer) to the integer heights. From we see that for half-odd- integer S 
the height is restricted to be an integer on one of the four sublattices, an in- 
teger plus 1/4 on the second, an integer plus 1/2 on the third, and an integer 
plus 3/4 on the fourth; the fractional parts of these heights are as shown in 
Fig&; the steps between neighboring heights are always an integer plus 1/4, 
or an integer plus 3/4. For S an odd integer, the heights are integers on one 
square sublattice, and half-odd-integers on the second sublattice. Finally for 
even integer S the offset has no effect and the height is an integer on all sites. 
We discuss these classes of S values in turn in the following subsections. 

4.1.1 S even integer 

In this case the offsets 2SXj are all integers, and 150(1 is just an ordinary 
three dimensional height model which has been much studied in the literature 
| 57l I59| . Unlike the two-dimensional case, three-dimensional height models 
generically have no roughening transition, and the interface is always smooth 
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\/ /\ \/ /\ 
/\ V /\ \/ 

w-^x^w-^x 

\/ /\ \/ A 

Z-^Y^Z-^'' 
(c) 

Fig. 9. Specification of the non-zero values of the fixed fields (a) Xj, (b) J^j^, (c) 
t^vxAvVix introduced in 14911 . The notational conventions are as in Fig |H| Only 
the fi — T components of yj^j, are non-zero, and these are shown in (b). Only the 
spatial components of efivxA^yjx are non-zero, and these are oriented as in (c) with 
magnitude 1/4. The four dual sublattices, W, X, Y, Z, are also indicated in (c). 
Note that Xw = 0, Xx = 1/4, Xy = 1/2, and Xz = 3/4. 



| 59) . With all heights integers, the smooth phase breaks no lattice symmetries. 
So square lattice antiferromagnets with S even integer can have a paramag- 
netic ground state with a spin gap and no broken symmetries. The smooth 
interface corresponds to confinement in the dual compact U(l) gauge theory 
|6U|: consequently the Za of Z are confined, and the elementary excitations 
are S = 1 quasiparticles, similar to the ipa of Sip. This is in accord with the 
exact ground state for a 5 = 2 antiferromagnet on the square lattice found 
by Affleck et al, the AKLT state EI]. 



4.1.2 S half-odd-integer 



Now the heights of the interface model can take four possible values, which 
are integers plus the offsets on the four square sublattices shown in Fig|5Ji. 
As in Section ETI l. the interface is always smooth i.e. any state of H5U|I has a 
fixed average interface height 
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^-]v:E(^.)' (52) 

where the sum is over a large set of Nd dual lattice points which respect the 
square lattice symmetry. Any well-defined value for H breaks the uniform 
shift symmetry of the height model under which Hj — )■ Hj± 1. In the present 
context, only the value of H modulo integers is physically significant, and so 
the breaking of the shift symmetry is not important by itself. However, after 
accounting for the height offsets, we now prove that any smooth interface must 
also break a lattice symmetry with the development of bond order, this means 
that Za in (|44(l describes spin gap ground states of the lattice antiferromagnet 
which necessarily have spontaneous bond order. 

The proof of this central result becomes clear upon a careful study of the 
manner in which the height model in (|5()|l and H51|l implements the 90° rotation 
symmetry about a direct square lattice point. Consider such a rotation under 
which the dual sublattice points in FiglSJ: interchange as 

W ^X, X ^Y, Y ^ Z, Z ^W. (53) 

The terms in the action in 151() will undergo a 90° rotation under this trans- 
formation provided the integer heights hj transform as 

hw hx, hx hy, hy hz, hz ^ hw ~ 1- (54) 

Notice the all important —1 in the last term — this compensates for the 'branch 
cut' in the values of the offsets Xj as one goes around a plaquette in FigO;. 
From it is evident that the average height H H — 1/4 under the 90° 
rotation symmetry under consideration here. Hence, a smooth interface with 
a well-defined value of H always breaks this symmetry. 

We now make this somewhat abstract discussion more physical by pre- 
senting a simple interpretation of the interface model in the language of the 
S = 1/2 antiferromagnet 62 . From Fig it is clear that nearest neigh- 
bor heights can differ either by 1/4 or 3/4 (modulo integers). To minimize 
the action in (|50|l . we should choose the interface with the largest possible 
number of steps of ±1/4. However, the interface is frustrated, and it is not 
possible to make all steps ±1/4 and at least a quarter of the steps must be 
±3/4. Indeed, there is a precise one-to-one mapping between interfaces with 
the minimal number of ±3/4 steps (we regard interfaces differing by a uni- 
form integer shift in all heights as equivalent) and the dimer coverings of the 
square lattice: the proof of this claim is illustrated in Fig^| We identify each 
dimer with a singlet valence bond between the spins (the ellipses in Figl^J, 
and so each interface corresponds to a quantum state with each spin locked 
in a singlet valence bond with a particular nearest neighbor. Fluctuations of 
the interface in imaginary time between such configurations correspond to 
quantum tunneling events between such dimer states, and an effective Hamil- 
tonian for this is provided by the quantum dimer model [631 . While such an 
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Fig. 10. Mapping between the quantum dimer model and the interface model 2^ 
in I50I I. Each dimer on the direct lattice is associated with a step in height of ±3/4 
on the link of the dual lattice that crosses it. All other height steps are ±1/4. Each 
dimer represents a singlet valence bond between the sites, as in Fig|5] 

interpretation in terms of the dimer model is appealing, we should also note 
that it is not as general as the dual interface model: on certain lattices, while 
the collinear paramagnetic state continues to have a representation as a dual 
interface model, there is no corresponding dimer interpretation |64| . 

The nature of the possible smooth phases of the interface model are easy to 
determine from the above picture and by standard techniques from statistical 
theory |2l)[ I62| . As a simple example, the above mapping between interface 
heights and dimer coverings allows one to deduce that interfaces with average 
height H = 1/8,3/8,5/8,7/8 (modulo integers) correspond to the four-fold 
degenerate bond-ordered states in Fig lllb . To see this, select the interface 
with hj = for all j: this interface has the same symmetry as Figlllb. and 
a simple computation summing over sites from (|51|l shows that this state has 
average height H = -(0 +_l/4 + 1/2 + 3/4)/4 = -3/8 for S = 1/2. The 
remaining three values of H correspond to the three other states obtained 
by successive 90° rotations of Fig lllh ,. In a similar manner, interfaces with 
H = 0, 1/4, 1/2,3/4 (modulo integers) correspond to the four-fold degenerate 
plaquette bond-ordered states in Fie lTTb. A simple example of such an inter- 
face is the "disordered-flat" state ||65| in which /ij = on all sites j, except for 
the W sublattice which have Xj = 0; for these sites we have hj fluctuate ran- 
domly between hj = and hj = 1, and independently for different j. The aver- 
age height of such an interface is!H'=-((0 + l)/2 + l/4 + l/2 + 3/4)/4=-l/2 
for S = 1/2, and the mapping to dimer coverings in Fig 1101 shows easily 
that such an interface corresponds to the state in Fig lllb. All values of H 
other than those quoted above are associated with eight-fold degenerate bond- 
ordered states with a superposition of the orders in Fig lllb and b. 

All these phases are expected to support non-zero spin quasiparticle exci- 
tations which carry spin S = 1, but not S = 1/2. Despite the local corrugation 
in the interface configuration introduced by the offsets, the interface remains 
smooth on the average, and this continues to correspond to confinement in the 
dual compact U(l) gauge theory [Snj. Consequently the spinous of Fig|3|3 are 
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(a) (b) 

Fig. 11. Sketch of the two simplest possible states with bond order for 5* = 1/2 
on the square lattice: (a) the columnar spin-Peierls states, and (b) plaquette state. 
Here the distinct line styles encode the different values of the bond order parameter 
Qij in 12911 on the links. This should be contrasted from Figs. [Till where the line 
styles represented distinct values of the exchange constants in the Hamiltonian. In 
the present section, the Hamiltonian has the full symmetry of the square lattice, 
and the orderings represented above amount to a spontaneous breaking of the lat- 
tice symmetry. Both states above are 4-fold degenerate; an 8-fold degenerate state, 
with superposition of the above orders, also appears as a possible ground state of 
the generalized interface model. Numerical studies of a number of two-dimensional 
quantum antiferromagnets [H^ |571 EHl El ESI EI ESI have found ground states with 
spontaneous bond order, similar to the states shown above. 



confined in pairs. The structure of the resulting 5=1 triplon quasiparticles 
is very similar to the excitations of the paramagnetic phase of the coupled 
dimer antifcrromagnet of Section |21 as we already noted in Section ^ 

Support for the class of bond-ordered states described above has appeared 
in a number of numerical studies of S" = 1/2 antiferromagnets in d = 2 
which have succeeded in moving from the small g Neel phase to the large 
g paramagnet. These include studies on the honeycomb lattice |66| (duality 
mapping on the honyecomb lattice appears in Ref. ^19i), on the planar py- 
rochlore lattice |^ EHI (duality mapping for a lattice with the symmetry of 
the planar pyrochore is in Refs. [691 154| . with a prediction for the bond order 
observed), on square lattice models with ring-exchange and easy-plane spin 
symmetry '70' (duality mapping on spin models with easy plane symmetry is 
in Refs. 71, 52, 72 ), and square lattice models with SU(Af) symmetry EHl 
(the theories (|40|l . with a = \ . . . N , and H5U|I apply unchanged to SU(A'') an- 
tiferromagnets). The case of the square lattice antifcrromagnet with first and 
second neighbor exchange is not conclusively settled: while two recent studies 
I^Cl] (and earher work j5Sl|7B]) do observe bond order in a paramagnetic 
spin-gap state, a third has so far not found such order. It is possible that 
this last study is observing signatures of the critical point between the Neel 
and bond-ordered states (to be described in Section ^3 which is expressed 
in a theory for deconfined spinous in Zc in H55|l . 
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Finally, we also mention that evidence for the spontaneous bond order of 
Fig llll appears in recent numerical studies of doped antiferromagnets |78[I79 |. 

4.1.3 S odd integer 

This case is similar to that S half-odd-integer, and we will not consider it in 
detail. The Berry phases again induce bond order in the spin gap state, but 
this order need only lead to a two-fold degeneracy. 

4.2 Critical theory 

We turn finally to the very difficult issue of the nature of the quantum phase 
transition from the Neel state to one of the bond-ordered paramagnetic states 
in Fig IIUI as a function of increasing g. This has been a long-standing open 
problem, and many different proposals have been made. The two phases break 
different symmetries of the Hamiltonian, and so are characterized by very 
different order parameters (one lives in spin space, and the other in real space). 
Landau- Ginzburg- Wilson (LGW) theory would imply that a generic second- 
order transition is not possible between such phases, and one obtains either 
a first-order transition or a region of co-existence of the two orders. However, 
the bond-order in the paramagnet was obtained entirely from quantum Berry 
phases attached to the fluctuating Neel order, and it is not clear that LGW 
theory applies in such a situation. 

Recent work by Senthil et al. jSOl 1^ has proposed an elegant resolution 
to many of these problems, and we will describe their results in the remainder 
of this subsection. The results are based upon solutions of a series of simpler 
models which strongly suggest that related results also apply to the SU(2) 
invariant, S = 1/2 models of interest. The computations are intricate, but the 
final results are quite easy to state, and are presented below. We will mainly 
limit our discussion here to the case of antiferromagnets of spin S = 1/2. 

First, contrary to the predictions of LGW theory, a generic second-order 
transition between the Neel state and the bond-ordered paramagnet is indeed 
possible (let us assume it occurs at g = gc for Z in H4()|l l. The theory for such a 
quantum critical point is obtained simply by taking a naive continuum limit of 
Z while ignoring both the compactness of the gauge field and the Berry phases. 
Remarkably, these complications of the lattice model Z, which we have so far 
stated were essential for the complete theory, have effects which cancel each 
other out, but only at the critical point. Note compactness on its own is a 
relevant perturbation which cannot be ignored i.e. without Berry phases, the 
compact and non-compact lattice CP^ model have distinct critical theories 
1221. However, the surprising new point noted by Senthil et al. [501 is that 
the non- compact CP^ model has the same critical theory as the compact CP^ 
model with S = 1/2 Berry phases. Taking the naive continuum limit of Z in 
ipnjl . and softening the hard-constraint on the Zja, we obtain the proposed 
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theory for the quantum critical point between the Neel state and the bond- 
ordered paramagnet for spin S — 1 /2 [ 



-Zc = y" Vzair, t)VA^ (r, r) exp | - j d^rdr 



\{d^^lA^)Za\' + s\Zaf 



If 



(55) 



We have also included here a kinetic term for the A^, and one can imagine 
that this is generated by integrating out large momentum Zja- On its own, 
describes the transition from a magnetically ordered phase with Za condensed 
at s < Sc, to a disordered state with a gapless U(l) photon at s > Sc (here 
Sc is the critical point of Zc). Clearly the s < Sc phase corresponds to the 
Neel phase of Z in (|40|l for g < gc- However, the s > Sc phase does not 
obviously correspond to the g > gc bond-ordered, fully gapped, paramagnet 
of Z. This is repaired by accounting for the compactness of the gauge field 
and the Berry phases: it is no longer possible to neglect them, while it was 
safe to do so precisely at g = gc- The combined effects of compactness and 
Berry phases are therefore dangerously irrelevant at g — ge- 
lt is important to note that the critical theory of H55|) is distinct from 
the critical theory S^p in H12|l. although both theories have a global 0(3) 
symmetry |82j . In particular the values of the exponents v are different in the 
two theories, and the scaling dimension of the Neel order parameter (pa under 
Sip is distinct from the scaling dimension of the Neel order parameter z'^a'^^Zi, 
at the critical point of Zc- 

It is interesting that Zc in H55|) is a theory for the S* = 1/2 spinors Za- These 
can be understood to be the continuum realization of the spinous shown earlier 
in Fig|2|3. Thus the spinons become the proper elementary degrees of freedom, 
but only at the quantum critical point. Hence it is appropriate to label this 
as a 'deconfined quantum critical point' (80) . These spinons are confined into 
a 5* = 1 quasiparticle once bond order appears for g > gc, for reasons similar 
to those illustrated in Figl^J). 

A key characteristic of this 'deconfined' critical point is the irrelevance of 
the compactness of the gauge field, and hence of monopole tunnelling events. 
A consequence of this is that the flux of the gauge field in Zc is conserved. 
This emergent conservation law, and the associated long-range gauge forces 
are key characteristics of such critical points. 

We summarize in Fig 1121 our results for S* = 1/2 square lattice antiferro- 
magnets, as described by Z in (|40|l . 

The claims above for the conspiracy between the compactness and Berry 
phases at the critical point are surprising and new. They are central to a com- 
plete understanding of square lattice antiferromagnets, and a full justification 
of the claims appears in the work of Senthil et al. The following subsections 
illustrate their origin by considering a series of models, of increasing complex- 
ity, where similar phenomena can be shown to occur. The reader may also 
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Fig. 12. Phase diagram of the model Z in II4UI I of S = 1/2 antiferromagnets with 
full square lattice symmetry. There is a Neel phase for g < Qc which breaks spin 
rotation invariance; it has a doublet of gapless spin wave excitations. The bond- 
ordered paramagnet for g > Qc preserves spin rotation invariance but breaks square 
lattice symmetry; it has a gap to all excitations, and the non-zero spin excitations 
are described by S = 1 triplet quasiparticles which are very similar to the 'triplons' 
discussed in Section [2.11 The critical point at g = gc is described by the theory of 
S — 1/2 'spinous', Zc in 15511 at its critical point s — Sc; note that this mapping 
to the spinon theory Zc does not work away from g ~ gc, and spinous are confined 
for all g > gc- A phase diagram like the one above has been used as a point of 
departure to obtain a phase diagram for doped Mott insulators |22l I83| . as a de- 
scription of the cuprate superconductors; evidence for spontaneous bond order in 
doped antiferromagnets appears in Refs. |78II79| . 



find it useful to look ahead to Tables 1 and 2, which summarize the intricate 
relationships between the models considered. 

4.2.1 Lattice model at AT = 1 

This subsection describes a simplified lattice gauge theory model introduced 
by Sachdev and Jalabert While the duality analysis presented below was 
initiated in Ref. |51| . its correct physical interpretation, and the implications 
for more general models are due to Senthil et al. |8(J[ I5T] . 

The model of interest in this subsection is the = 1 case of Z. Physically, 
such a model will be appropriate for an antiferromagnet in the presence of a 
staggered magnetic field: such a field will prefer Z| over (say). So we write 
the preferred single component complex scalar simply as Zj = e'^^ , and obtain 
from ignj) 
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2i - n S ^ -p (^1 X: cos ie,.,A.A,,) 

3,t^ 3 / 



(56) 



We have chosen here to exphcitly include a compact Maxwell term for the 
gauge field, as that proves convenient in the description of the duality map- 
pings. Note that if we integrate out the 9j for large g, then we again obtain 
the model Za in l|43|) which was used to describe the paramagnetic phase in 
Section E!n So bond order appears also in the model Zi at large g. This bond 
order disappears as g is reduced, at a transition we will describe below. 

Rather than attack Zi directly, it is useful as a warm-up, and to make 
contact with previous work, to consider a sequence of simpler models that 
have been considered in the literature. As we have emphasized, Zi features 
the combined complications of compactness and Berry phases, essential for 
a proper description of quantum antiferromagnets. It is the simplest model 
in which it can be shown that these complications effectively neutralize one 
another at the critical point. 

In the following subsection, we make things simpler for ourselves momen- 
tarily by dropping both the compactness and the Berry phases. We will then, 
in the subsequent subsections, add these complications back in. 

A. XY model with a non-compact U(l) gauge field 

Dropping both compactness and Berry phases, Zi reduces to 

+ -Y,cos{A^,9j~A,M 



(57) 



Notice that the Maxwell term for the gauge field now has a simple Gaus- 
sian form. This is simply the lattice, classical, Ginzburg-Landau model (or 
an XY model) of a superconductor at finite temperatures coupled to electro- 
magnetism. This model has been studied extensively in the past, and the key 
result was provided by Dasgupta and Halperin j84| . As we review below, they 
showed that Za,c exhibited an inverted XY transition i.e. it was dual to the 
theory of a complex scalar ip in the absence of a gauge field: 

2sc,duai = y"pi^(r,r)exp (^-y"dVdr (^|9^,/-p + s|i/^|2 + :^|V7|4^^ (58) 

The field ^ is a creation operator for vortices in the original theory of the 
Ginzburg-Landau superconductor. These have a short-range interaction {u 
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above) because of the screening provided by the electromagnetic flux quantum 
attached to every vortex in (|57|) . So the vortex loops of (|57|l behave like the 
world lines of the dual boson field of H58|) . The tuning parameter s in H58|l 
is 'inverted' from the perspective of the direct theory: the s < Sc phase with 
(tp) 7^ has a vortex condensate and so is the normal state of a Ginzburg- 
Landau superconductor, while the s > Sc phase with {ip) = has the vortices 
gapped as in the superconducting phase. 

We now provide a few steps in the analysis which links H57|l to (|58|l . The 
steps are very similar to those described in Section 1^1 below (|43ll and (|44|l . 
We write the cosine in H57|l in its Villain form, decouple it by the Poisson 
summation formula using integer currents Jj^ , and also decouple the Maxwell 
term by a Hubbard-Stratonovich field Pj^; this yields the analog of (|45|l for 



J/J 



The advantage of this form is that the integrals over 9j and Aj^^ can be 
performed exactly, and they lead to the constraints 

^f^J]f^ = ; Jjf^ = e^^xAiyPjx. (60) 
We solve these constraints by writing 

Jjfi = e^^A^^&jA ; Pjfi = bj^ - ^iiV], (61) 

where bj^ is an integer valued field on the links of the dual lattice, and ipj is 
a real valued field on the sites of the dual lattice. This transforms (|59|) to 

Zsc,2 = n / ^"^J E "y E ^ ^t^flf 

-|5](e^.,Zi,6,A)'j; (62) 

precisely this dual form was obtained by Dasgupta and Halperin [SJ , and used 
by them for numerical simulations. We proceed further analytically, using 
methods familiar in the theory of duality mappings [S7j: we promote the 
integer valued to a real field by the Poisson summation method, and 
introduce, by hand, a vortex fugacity j/t,. This transforms Zsc,2 to 



/C30 f^OO f^OQ / 2 

%M / dipj / diJjexp {bjf, - 

-OO 'J —OO OO \ J fi 
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(63) 




Notice that the effect of the vortex fugacity is to yield the least action when 
bj^ is an integer (ignore i?j momentarily): so we have effectively 'softened' 
the integer constraint on bj^. We have also introduced here a new real valued 
field on the sites of the dual lattice simply to make the Zsc,3 invariant 
under U(l) gauge transformations of bj^. This is mainly because the physics 
is clearer in this explicitly gauge-invariant form. We could, if we had wished, 
also chosen a gauge in which = 0, and then the field would not be present 
in Zsc,3 (this justifies neglect of above). In the complete form in H63|l . it 
is clear from the first two Gaussian terms that fluctuations of the fej^ gauge 
field have been 'Higgsed' by the real field ipj. Indeed, it is more convenient to 
choose a gauge in which ipj = 0, and we do so. Now the fluctuations of 6jp 
are 'massive' and so can be safely integrated out. To leading order in z/y, this 
involves simply replacing with the saddle point value obtained from the 
first two Gaussian terms, which is 6jp — 0. So we have the very simple final 
theory 



which has the form of the dual XY model. We now take the continuum limit of 
l|64(l by a standard procedure (HEj of introducing a complex field ijj conjugate to 
e"'^', and obtain the theory Zgc,duai as promised. This establishes the duality 
mapping of Dasgupta and Halperin 84 . 

B. XY model with a compact U(l) gauge field 

Now we ease towards our aim of a duality analysis of Zi , by adding one layer 
of complexity to Za,c- We make the gauge field in l|57|l compact by including 
a cosine Maxwell term 



The Dasgupta-Halperin duality mapping can be easily extended to this theory. 
We now write both cosine terms in their Villain forms, and then proceed as 
described above. The results (|59|l and H62() continue to have the same form, 
with the only change being that the fields Pj^ and Lpj are now also integer 
valued (and so must be summed over) . Promoting these integer valued fields to 
real fields by the Poisson summation method following Ref. we now have 




(64) 
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to introduce two fugacities: a vortex fugacity ?/„ (as before), and a nionopole 
fugacity (discussed below). Consequently, Zgc,3 in (|63|l now takes the form 



/oo poo i^oa / g2 

dbj^ / d^j / diJjexp -y XI " ^m'^j)^ 
J -oo J~co J-oc \ ^ 

- I X {(■^ux^uhjxf' +yv'Y COS (27r6jp - A^^j) 



+yrn Y COS {2n(pj - l?j)^ 



(66) 



Again, the positions of the iJj above are dictated by gauge invariance, and 
the effect of the vortex and monopole fugacities is to soften the integer value 
constraints on the and ipj. Proceeding as described below (|63|) . we work 
in the gauge (pj — 0, and to leading order in yy, y„i replace by its saddle 
point value in the Gaussian part of the action, which remains = 0. Then, 
instead of l|64|l . we obtain 



Zma — 



Yl / dt^jexp X cos (A^dj) + y^ ^ cos (i?j) . (67) 



We see that the new second term in (|67|l acts like an ordering field on the dual 
XY model. Taking the continuum limit as was done below H64|l using BS; a 
complex field ^ conjugate to e*''^, now instead of -Zsc.duai in we obtain 

IH7IIHH1 

ZM,duai = j I?^(r,T)exp|^- j d2rdr(^|5,,Vp + s|7/;|2 

+ f|V'|'-2/m(^ + V*))^ (68) 

The new term proportional to y„i has the interpretation of a monopole fugac- 
ity. The compact gauge field now permits Dirac monopoles, which are points 
in spacetime at which vortex loops of the 'superconductor' can end: hence 
is coupled to the creation and annihilation operators for the dual boson ip i.e. 
the vortices. In the form H68|l it is also clear that ym acts like an ordering 
field in the dual XY model. We expect that such an XY model has no phase 
transition, and {tp) ^ for all s. So the presence of monopoles has destroyed 
the 'superconducting' phase. Comparing the properties of (|58|l and H68() we 
therefore conclude that making the gauge field compact in Z^c in (|57|l is a 
strongly relevant perturbation: the inverted XY transition of Z^c is destroyed 
in the resulting model Zm- 
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C. Berry phases 

We are finally ready to face Zi, and add in the final layer of complication of 
the Berry phases. Again, the Dasgupta-Halpcrin duality can be extended by 
combining it with the methods of Section 14.11 (this was partly discussed in 
Ref. O]). Now the monopoles carry Berry phases |89[ I19| . and these lead to 
cancellations among many monopole configurations. In the long- wavelength 
limit it turns out that the only important configurations are those in which 
the total monopole magnetic charge is q times the charge of the elementary 
monopole ISHlEll^ni- Here q is the smallest positive integer such that 



i.e. q — A for S half an odd integer, q = 2 for S an odd integer, and q — 1 
for S an even integer. Using the physical interpretation of (|68|l . we therefore 
conclude that the monopole fugacity term should be replaced by one in which 
the monopoles are created and annihilated in multiples of q; the dual theory 
of Zi in (|56|l then becomes 



An explicit derivation of the mapping from Zi to Zi^^uai can be obtained 
by an extension of the methods described above for Zsc and Zm ■ We express 
the Berry phase term using the 'background field' a?^ in 146(1 . and then we 
find that 2^sc,2 in H62|l is now replaced by 



Notice that, as in Section 1^1 the Berry phases appear as offsets in the dual 
action. We now promote the integer field and cpj to real fields by the 
Poisson summation method (just as in JHSJ), at the cost of introducing vortex 
and monopole fugacities. The final steps, following the procedure below (|66|l . 
are to transform to the gauge cpj = 0, and to then set the 'Higgsed' dual 
gauge field to its saddle point value determined from the Gaussian terms 
in the action. It is the latter step which is now different, and the presence 
of the now implies that the saddle point value bj^ will be non-zero and 
site dependent. Indeed, it is crucial that the saddle point be determined with 
great care, and that the square lattice symmetry of the underlying problem be 



(69) 
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fully respected. This saddle point determination is in many ways analogous 
to the computation in Section III.B of Ref. [20], and it is important that 
all the modes on the lattice scale be fully identified in a similar manner. 
The similarity to Ref. |2Dj becomes clear after using the parameterization in 
H49() for aj^ in terms of the Xj and the 3^j^ shown in Fig |2| Finally, after 
transforming hj^ + 2SA^Xj and ??j -dj + ATTSXj, we obtain from 1)7111 



Now, the saddle point value of the massive field bj^ is easily determined from 
the first terms in H72|l. yielding 



where a = 25'e^/(e^ + 8g). Note that only the spatial components of are 
non-zero, and these have the simple structure of Fig In particular, the 
magnitude of the bj^ are the same on all the spatial links, and the use of 
was crucial in obtaining this appealing result. With this saddle point value, 
H72|) simplifies to the following model for the field only (this is the form of 
Ht)7|) after accounting for Berry phases): 



The most important property of this dual XY model is the nature of the 
ordering field in the last term of (|74|) . For S = 1/2, notice from Fig|5^ that 
this field is oriented north/east /south/ west on the four sublattices in of the 
dual lattice in Fig So if we take a naive continuum limit, the average 
field vanishes! This is the key effect responsible for the cancellations among 
monopole configurations induced by Berry phases noted earlier; in the dual 
formulation, the Berry phases have appeared in differing orientations of the 
dual ordering field. The XY model in (|74(l also has the contribution from 
which appear as a 'staggered flux' acting on the iJj (see FiglHl:;), but we now 
show that this is not as crucial in the continuum limit. 

Before we take the continuum limit of -Zi,4, we discuss its implementa- 
tion of the square lattice symmetries. In particular, we are interested in the 




(72) 



(73) 
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^4 symmetry which rotates the four sublattices in Fig |^ into each other, 
as the values of Xj seem to distinguish between them. Let us consider the 
symmetry 7?.„ which rotates lattice anticlockwise by an angle mr/2 about the 
direct lattice point at the center of a plaquette in FigO;, associated with the 
transformation in H53() . It is easy to see that Zi^4 remains invariant under TZn 
provided we simultaneously rotate the angular variables i!)f. 

Tin ■■ i)j + uSt: (75) 

It is now useful to introduce complex variables which realize irreducible rep- 
resentations of this Z4 symmetry. We divide the lattice into plaquettes like 
those in Fig|5t, and for each plaquette we define variables ipp, with p integer, 

by 

Note that we need only usep = 0, 1, 2, 3 because ipp depends only onp(mod 4). 
Under the symmetry TZ^ we clearly have 

Tin-. Vp ^ e™(2^-P)"/Vp; (77) 

the factor of e*"'^'^ arises from (|75|l . and that of from the real-space 

rotation of the lattice points. Note that only for p = 25 is ipp invariant under 
TZn, and this is consistent with the fact that it is ip2S which appears in Zi_4 as 
the ordering field term. Let us now write the action in 2i 4 in terms of these 
new variables. Ignoring the spacetime variation from plaquette to plaquette, 
the action per plaquette is 

3 

Si,4 = -22/„ ^ [cos {tt{p - a)/2) iV'pl'] - Vm ii^2S + i^h) + ■■■ (78) 

p=0 

Here the ellipses represent other allowed terms, all consistent with the sym- 
metry ()77|l . which must be included to implement the (softened) constraints 
on i/jp arising from H76() and the fact that the e^^^ are uniniodular. Apart from 
'025, for which there is already an ordering field in the action, the condensa- 
tion of any of the other ipp breaks the lattice symmetry (I77II . and so drives a 
quantum phase transition to the bond-ordered state. The choice among the ipp 
is controlled by the co-efficient of the ?/„ term in (|78|l . and we choose the value 
oi p ^ 2S for which cos (7r(a + p)/2) is a maximum. We are interested in the 
large g paramagnetic phase, and here a is small, and the appropriate value 
is p = 0. The resulting continuum theory for ip = ipo then must be invariant 
under (|77|l . and it is easily seen that this has just the form Zi.duai in H7U|) 
with q determined by (|69|l . Other choices of p for the order parameter lead to 
different types of bond order, with a ground state degeneracy smaller or larger 
than the q in (|69|l : such states have partial or additional bond order, and are 
clearly possible in general. However, our analysis of the paramagnetic states 
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iV = 1 



Direct lattice model 


Dual model 


Csc = (1/(26^)) (e^uxA,A,xf 


•CsCdual 




>' 


Cm = -(1/e^) cos (e^^A/ii/^jA) 
- (l/g)cos(Z\^6l, 


•^^Af.dual 






Ci = -(l/e^)cos(e^„A'4„AjA) 

- (1/g) cos {A^e, - A,^) - i2Sr,jA,r 


-Cl,dual 







Table 1. Summary of the duality mappings for TV = 1. Only the Lagrangean's 
are specified, and a summation/integration of these over spacetimes is implicit. The 
fixed field rjj = ±1 in the Berry phase in the third row is the sublattice staggering 
factor in |(2J. The integer q in the third row is specified in For S = 1/2, we have 
g = 4, and then the ymq perturbation is dangerously irrelevant. Hence the critical 
theory for the model with monopoles and Berry phases in the third row, is identical 
to that for the first row. 

in Section 14.11 indicates that a choice = V'pt^o is unhkely for the models 
under consideration here, and we will not consider this case further here. 

We have now completed our promised derivation of the model 2^i,duai in 
l|7()(l dual to the = 1 lattice gauge theory model Zi in (jHEJ. Rather than 
being an XY model in a field (as in Hfi8|l 'l. Zi^duai is an XY model with a q- 
fold anisotropy. This anisotropy encapsulates the g-fold binding of monopoles 
claimed earlier. In the language of (|74() the average ordering fields on the dj 
oscillate from site to site and cancel out, and only the q-ih moment of the field 
survives. Now the combined effect of the monopoles and Berry phases in Zi 
is decided by the term proportional to Hmq- In the paramagnetic phase of the 
direct model, which is s < Sc and (V') ^ 0, this q-fold anisotropy is certainly 
very important. For S — 1/2, q = 4 it orders the ip field along four particular 
angles, and these are easily shown to be [SJ one of the four degenerate bond- 
ordered states in Fig However, at the critical point s = Sc it is known 
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that this 4-fold anisotropy is irrelevant 90 : so in Zi juai the monopoles can 
be neglected at the critical point s = Sc, but not away from it. 

We have now achieved the desired objective of this subsection. Compact- 
ness alone was a strongly relevant perturbation on the model of a scalar field 
coupled to electromagnetism in Zsc- However, when we combined compact- 
ness with the Berry phases in Zi , then we found that the monopoles effectively 
cancelled each other out at the critical point for S = 1/2. Consequently the 
theory for the critical point in Zi is identical to the theory for the critical 
point in Zgc, and this is the simple inverted XY model 2^sc,duai in (|58|l . The 
results of this subsection are summarized in Table 1. 

4.2.2 Easy plane model at TV = 2 

A second explicit example of the remarkable phenomenon described above 
is provided by the physically relevant N = 2 case of the model of central 
interest, Z in H40|l . but in the presence of an additional spin-anisotropy term 
preferring that the spins lie within the XY plane. In such a situation, we may 
write the complex spinor Zja as 

"^" = 7lU^'^v' ^^^^ 

so that the action is expressed in terms of two angular fields, 0| and O^. 
Inserting H79|l in (|40|) , we obtain a generalization of the = 1 model Zi in 

+ Y i^^t^^ja - Ajf,) + i2SYv3Ajr ■ (80) 

As in (|56|l . we have chosen to explicitly include a Maxwell term for the U(l) 
gauge field as it proves convenient in the subsequent duality analysis. The 
model Z-i provides a complete description of the phases of the square lattice 
antiferromagnet (|28|l with an additional easy-plane anisotropy term. 

We can now proceed with a duality analysis of H80|) using methods precisely 
analogous to those discussed in Section [4.21 1: the only difference is we now 
have two angular fields Qa=\,i^ and so certain fields come with two copies. We 
will therefore not present any details, and simply state the series of results 
which appear here, which closely parallel those obtained above for A^ = 1. 

• Neglecting both compactness of the C/(l) gauge field and the Berry phases, 
it is straightforward to take the continuum limit of Z2 in its direct rep- 
resentation, and we obtain the theory Z^ in IjSSI) . but with an additional 
spin-anisotropy term 
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J Vza{r,T)VA^{r,T)exp(^^ J dPrdr 



(81) 



where u > prefers spins in the easy plane. We can carry through the 
analog of the duality mapping between H57|l and H58|l . and instead of H58|) 
we now obtain a theory for two dual fields ipa representing vortices in 9^ 
and 6*1 [12 



2c,dual 



j VMr, T)VB^{r, r) exp f - ^ d'^rdr 



\{d,-iBM' 



+ \id^ + iBM^ + s\iba\^ + -(|V;,|2)2 +TJ|^^|2|V,J2 



(82) 



Note that there is now a non-compact U(l) gauge field which survives 
the continuum limit: this field arises from the analog of the field bj^ in 
and here it is not completely Higgsed out. The most remarkable property 
of H82() is that it is identical in structure to (|81|l : the actions are identical 
under the mapping Z| ^ , Z| ^ V'l: and B^. In other words, the 

theory Z2c is self-dual 82 . 

As in Section r4.2l l. we next make the gauge field compact, but continue 
to ignore Berry phases i.e. we perform a duality analysis on (|80|) . in the 
absence of the last term in the action. Now, instead of H68I) . (|82|l is modified 
to 



-'2M,dual 



J VMr. T)VB^{r, r) exp ^- j d^rdr 



(83) 



The last term represents the influence of monopoles, and these now have 
the effect of turning a V't vortex into a ipl vortex |82[ I80L IST| . Again, as 
m 

(EHl, the Um term in 183(1 is clearly a strongly relevant perturbation 
to Z2c,duai in H82|l . It ties the phases of V't and V'i to each other, so that 
(|83|l is effectively the theory of a single complex scalar coupled to a non- 
compact U(l) gauge field B^. However, we have already considered such a 
theory in the direct representation in H57|l . We can now move from the dual 
representation in H83|l back to the direct representation, by the mapping 
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-'2M.dual 



dT 



between H57|l and H58|) . This leads to the conclusion, finally, that the theory 
(|83|l is dual to an ordinary XY model. In other words, the theory Z2 in H8U|) 
without its Berry phase term is an XY model. However, this is precisely 
the expected conclusion, and could have been easily reached without this 
elaborate series of duality mappings: just integrating over Aj^ for large 
yields an XY model in the angular field — which represents the 
orientation of the physical in-plane Neel order. 

Finally, let us look at the complete theory Z2 ■ An explicit duality mapping 
can be carried out, and as in (|70|l . the action (|83|l is replaced by |71l I52L 

EOllHIl 

J VMr, T)VB^{r, r) exp^- j d^r. 

+ ^(eM^A9.i3A)'-y™,((V^TV^i)''+(^i>T*)') (84) 

where the integer q was defined in (|69|l . The subsequent reasoning is the 
precise analog of that for iV = 1. For 5 = 1/2 and 9 = 4, the term pro- 
portional to ymq representing g-fold monopole is irrelevant at the critical 
point (but not away from it in the paramagnetic phase) . Consequently, the 
critical theory of H84|) reduces to (|82|l . So just as at iV = 1, the combined 
influence of monopoles and Berry phases is dangerously irrelevant at the 
critical point, and for the critical theory we can take a naive continuum 
limit of Z2 neglecting both the Berry phases and the compactness of the 
gauge field. 

We have now completed our discussion of the N — 2 easy plane model and 
established the existence of the same remarkable phenomenon found in Sec- 
tion ^21 1 for = 1, and claimed more generally 80, 81 at the beginning of 
Section as the justification for the critical theory H55|) . As we saw in some 
detail in Section Wa\ monopoles, and attendant Berry phases, are absolutely 
crucial in understanding the onset of confinement and bond order in the para- 
magnetic phase. However, for S — 1/2, the Berry phases induce a destructive 
quantum interference between the monopoles at the quantum critical point, 
leading to a critical theory with 'deconfined' spinous and a non-compact U(l) 
gauge field which does not allow monopoles. These results are summarized in 
Table 2. 

The results in Tables 1 and 2 can be generalized to arbitrary values of 
N, for models with the analog of an 'easy plane' anisotropy: as in (|79|l . all 
the Za have equal modulus and are expressed in terms of a = 1 ... A angles 
9a- The dual models have A vortex fields ipa, and A^ — 1 non-compact U(l) 
gauge fields Bh^, b = I . . .{N — 1). For a = 1 ... (A^ — 1), the field -^a has a 
charge -t-1 under the gauge field with b = a, and is neutral under all gauge 
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N = 2, easy plane 



Direct lattice model 


Dual model 


- {1/ g) cos (A^Oja ~ Ajf,) 


+ -^{<^t^i'\d^Bxf 


C2M = -{l/e'^)cos{ef^^xA„Ajx) 
- (l/(2ff))cos(A,e,„-A,„) 


'C2A/,dual = (9m - i-Bfi)l/'t 


jC2 = -(l/e^) COS {ef,t,x A,, Ajx) 
-{l/g) cos {A^Oja - Aj^) - i2Sr]jAjr 


'C2,dual = I (9m - ^Bfj,)lpff 

+ \(d^+iB,)^^f + s\i,af + ^{\t^aff 

-j/-,((V'T^i)'' + (^T>r)') 



Table 2. As in Table 1, but for the N — 2 easy plane case. The index a extends over 
the two values t, J,. Again for S — 1/2, g = 4, the critical theory for the third row 
is the same as that for the first row. The dual model in the second row is effectively 
the theory of a single complex scalar coupled to a non-compact U(l) gauge field BmI 
by the inverse of the duality mapping in the first row of Table f , this theory has a 
direct XY transition. 



fields with h ^ a. For a = N, the field i/jn, has a charge —1 under all N — 1 
gauge fields. (This gauge structure is similar to that found in 'moose' field 
theories [SJ-) The dual representation of the monopole operator is V'o, 
and this appears as the co-efficient of ym (notice that this operator is neutral 
under all the gauge fields). The q*'^ power of this operator appears as the 
coefficient of Umq- Note that the monopole operators involves a product of 
A'^ fields, and for large enough N, both ym and ymq can be expected to be 
irrelevant perturbations at the quantum critical point. 
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Finally, these analyses of Z in H4U|I can be complemented by a study of 
its — !■ oo limit, without any easy-plane anisotropy. This was carried out 
some time ago [921 151) , and it was found that monopoles were dangerously 
irrelevant at the quantum critical point, hoih with or without Berry phases 
(as noted above for large N in the easy plane case). It is important to note 
that the situation at large TV is subtly different from that for iV = 1, 2: in the 
latter case, monopoles are dangerously irrelevant in the presence of S" 1/2 
Berry phases, but relevant without Berry phases. The key understanding of 
this distinction emerged in the recent work of Senthil tt al. [SDl E| i which 
finally succeeded in placing the earlier large N results within the context of 
dual theories of topological defects in statistical mechanics. 

5 Triangular lattice antiferromagnet 

We continue our analysis of quantum antiferromagnets with an odd number 
of = 1/2 spins per unit cell, but consider a class qualitatively different from 
those in Section ^ One of the defining properties of the models of Section 0] 
was that the magnetically ordered Neel state was defined by the aver- 
age magnetic moment on all sites were collinear, and only a single vector n 
was required to specify the orientation of the ground state. This section shall 
consider models in which the moments are non- collinear, the triangular lat- 
tice is the canonical example. However, similar results should also apply to 
other two-dimensional lattices with non-coUinear ground states, such as the 
distorted triangular lattice found in CS2CUCI4 [TT] . 

We consider the model H28I) , but with the spins residing on the sites of the 
triangular lattice. This has a magnetically ordered state illustrated in Fig ll3l 
for this state ^ is replaced by 

(Sj) = No (ni cos(Q • r) + na sin(Q • r)) . (85) 

Here Q — 27r(l/3, 1/V3) is the ordering wavevector, and ni^2 are two arbi- 
trary orthogonal unit vectors in spin space 

nl = nl = l ; rii • ns = 0. (86) 

A distinct ground state, breaking spin rotation symmetry, is obtained for each 
choice of ni^2- 

We now wish to allow the values of ni^2 to fluctuate quantum mechanically 
across spacetime, ultimately producing a paramagnetic state. As in Section^ 
we should account for the Berry phases of each spin while setting up the effec- 
tive action: an approach for doing this is presented in Section VI of Ref . ■ 
However, the full structure of the critical theory is not understood in all cases, 
as we describe below. 

One possible structure of the paramagnetic state is a confining, bond- 
ordered state, similar to that found in Section^ However, there is no com- 
plete theory for a possible direct second-order transition from a non-coUinear 
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Fig. 13. Quantum phase transition described by in II9UI I as a function of s. The 
state on the left has non-collinear magnetic order described by II85|I . The state on the 
right is a 'resonating valence bond' (RVB) paramagnet with topological order and 
fractionalized S = 1/2 neutral spinon excitations (one spinon is shown above). Such 
a magnetically ordered state is observed in CS2CUCI3 |11II12) . and there is evidence 
that the higher energy spectrum can be characterized in terms of excitations of the 
RVB state El- 



magnetically ordered state to such a paramagnet. Ignoring Berry phases, one 
could define the complex field = niQ +m2a, which, by obeys = 0, 
and then proceed to write down an effective action with the structure of 
However, it is clear that such a theory describes a transition to a paramagnetic 
phase with a doublet of 5 = 1 triplet quasiparticles, and we can reasonably 
expect that such a phase has spontaneous bond order (in contrast to the ex- 
plicit dimerization in the models of Section Berry phases surely play an 
important role in inducing this bond order (as they did in Section l4.1|l . but 
there is no available theory for how this happens in the context of H14|l . Indeed, 
it is possible that there is no such direct transition between the non-collinear 
antiferromagnet and the bond-ordered paramagnet, and resolving this issue 
remains an important open question. 

In contrast, it is possible to write down a simple theory for a direct transi- 
tion between the non-collinear antiferromagnet and a paramagnetic phase not 
discussed so far: a resonating valence bond hquid with deconfined 

spinous and topological order. This theory is obtained by observing that the 
constraints (|86() can be solved by writing |94[ 



ni -I- in2 = eabWbO-acWc 



(87) 



where Wa is a 2 component complex spinor obeying -I- = 1. It is 

useful to compare H87(l with (|31|l . which parameterized a single vector also 
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in terms of a complex spinor Za- Whereas was invariant under the U(l) 
gauge transformation H32|) . notice that (|87|l is only invariant under the Z2 
gauge transformation 



Wair,T) — > g{r,T)wa{r,T) 



(88) 



where g{r,T) ~ ±1 is an arbitrary field which generates the gauge transfor- 
mation. This Z2 gauge transformation will play an important role in under- 
standing the structure of the paramagnetic phase 21.. .22. .96. .23^. .97j. 

We can now study fluctuations of the non-coUinear antiferromagnet by 
expressing the effective action in terms of the Wa- Apart from the familiar 
constraints of spin rotational invariance, and those imposed by (|88|l . the effec- 
tive action must also obey the consequences of translational invariance which 
follow from 1)85(1 : the action must be invariant under 



Wa WaB 



-iQ-a/2 



(89) 



where a is any triangular lattice vector. In the continuum limit, this leads to 
the following effective action 



d rdr 



(90) 



notice there is a free integration over the Wa, and so we have softened the rigid 
length constraint. Comparing this with (|55|l . we observe that the U(l) gauge 
field is now missing, and we simply have a Landau-Ginzburg theory for a 2 
component complex scalar. The Z2 gauge invariance H88|l plays no role in this 
continuum critical theory for the destruction of non-coUinear magnetic order, 
but as we discuss below, it will play an important role in the analysis of the 
paramagnetic phase. The theory ((9()|l has a global 0(4) invariance of rotations 
in the 4-dimensional space consisting of the real and imaginary parts of the 2 
components of Wa ■ consequently the critical exponents of 1)90(1 are identical to 
those of the well known 4-component (p^ field theory. Notice that there is no 
0(4) invariance in the microscopic theory, and this symmetry emerges only 
in the continuum limit 98, 95 : the simplest allowed term which breaks this 
0(4) invariance is \eabWadfi,Wb\'^ , and this term is easily seen to be irrelevant 
at the critical point of the theory ((90(1 . 

Let us now turn to a discussion of the nature of the paramagnetic phase 
obtained in the region of large positive s in ((90(1 . Here, the elementary exci- 
tations are free Wa quanta, and these are evidently S* = 1/2 spinous. There 
is also a neutral, spinless topological excitation 99, 2"l||22| whose importance 
was stressed in Ref. [HI]: this is the 'vison' which is intimately linked with 
the Z2 gauge symmetry ((88(1 . It is a point defect which carries Z2 gauge flux. 
The vison has an energy gap in the paramagnetic phase, and indeed across 
the transition to the magnetically ordered state. This was actually implicit 
in our taking the continuum limit to obtain the action l(90() . We assumed 
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that all important spin configurations could be described by a smooth, single- 
valued field Wa {r, t) , and this prohibits vison defects around which the Wa 
are double- valued. It is also believed that the vison gap allows neglect of 
Berry phase effects across the transition described by (|90f) : after duality, the 
Berry phases can be attached to monopoles and visons (23 I100| , and these 
are suppressed in both phases of FigllHI 

6 Conclusions 

This article has described a variety of quantum phases of antiferromagnetic 
Mott insulators, and the transitions between them. 

Let us first summarize the phases obtained in zero applied magnetic field, 
and transitions that can be tuned between them by varying the ratio of ex- 
change constants in the Hamiltonian (experimentally, this can be achieved by 
applied pressure). The magnetically ordered states discussed were the coUinear 
Neel state (shown in Figs El and II 2|l . and the non-coUinear 'spiral' (shown in 
Fig I13|l . We also foimd paramagnetic states which preserved spin rotation 
invariance and which had an energy gap to all excitations: these include the 
dimerized states (shown in Figs |2 and 0J, the related bond-ordered states 
which spontaneously break lattice symmetries (shown in Figs llll and ll2|l . and 
the 'resonating valence bond' paramagnet with topological order and decon- 
fined spinous (shown in Fig I13|l . The continuous quantum phase transitions 
we found between these states were: 

(a) the transition between the dimerized paramagnet and the coUinear Neel 
state (both states shown in FigEJ was described by the theory S^p in (|12|l : 
(6) the transition between the dimerized paramagnet and a non-coUinear mag- 
netically ordered state was described by in H14|) : 

(c) the transition between the coUinear Neel state and the paramagnet with 
spontaneous bond order (shown in Fig ll2(l was described for S = 1/2 antifer- 
romagnets by Zc in 1)55(1 : 

(d) the transition between the state with non-coUinear magnetic order and 
the RVB paramagnet (both states shown in Fig ll3|l was described by in 

We also mention here other quantum transitions of Mott insulators, which 
involve distinct paramagnets on both sides of the critical point. These we did 
not discuss in the present paper, but such transitions have been discussed in 
the literature: 

(e) the transition between a paramagnet with spontaneous bond order fFig ll2() 
and a RVB paramagnet fFig ll3|l is described by a compact U(l) lattice gauge 
theory with charge 2 Higgs fields (closely related to Zi in ((56|l l. and is dis- 
cussed in Refs. [TOnHTTm i^: 

(/) transitions between paramagnets with different types of spontaneous bond 
order can be mapped onto transitions between different smooth phases of 
height models like ((50|l . and are discussed in Refs. jl()2l 1103] . 
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Section 121 also considered quantum transitions that could be tuned by an 
applied magnetic field. We mainly considered the case of the coupled-dimer 
antiferromagnet, but very similar theories apply to the other states discussed 
above. The general theory has the structure of Sxp in H24|) . describing the Bose- 
Einstein condensation of the lowest non-zero spin quasiparticle excitation of 
the paramagnet. For the coupled dimer model this quasiparticle had 5 = 1, 
but an essentially identical theory would apply for cases with 5=1/2 spinon 
quasiparticles. 
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